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ABSTRACT 

The wavelength dependence of atmospheric refraction causes elongation of finite- 
bandwidth images along the elevation vector, which produces spurious signals in weak 
gravitational lensing shear measurements unless this atmospheric dispersion is cali- 
brated and removed to high precision. Because astrometric solutions and point spread 
function (PSF) characteristics are typically calibrated from stellar images, differences 
between the reference stars' spectra and the galaxies' spectra will leave residual er- 
rors in both the astrometric positions (Ai?) and in the second moment (width) of the 
wavelength-averaged PSF (AV) for galaxies. We estimate the level of AV that will 
induce spurious weak lensing signals in PSF-corrected galaxy shapes that exceed the 
statistical errors of the Dark Energy Survey (DES) and the Large Synoptic Survey Tele- 
scope (LSST) cosmic-shear experiments. We also estimate the AR signals that will 
produce unacceptable spurious distortions after stacking of exposures taken at different 
airmasses and hour angles. Using standard galaxy and stellar spectral templates we 
calculate the resultant errors in the griz bands, and find that atmospheric dispersion 
shear systematics, left uncorrected, are up to 6 and 2 times larger in g and r bands, 
respectively, than the thresholds at which they become significant contributors to the 
DES error budget, but can be safely ignored in i and z bands. For the stricter LSST re- 
quirements, the factors are about 30, 10, and 3 in g, r, and i bands, respectively. These 
shear systematic errors scale with observed zenith angle z as (tan 2 z), for which we 
take a nominal value of unity — simulations of DES and LSST suggest 0.6-1.0. We find 
that a simple correction linear in galaxy color is accurate enough to reduce dispersion 
shear systematics to insignificant levels in the r band for DES and i band for LSST, 
but still as much 5x above the threshold of significance for LSST r-band observations. 
More complex approaches to correction of the atmospheric dispersion signal will likely 
be able to reduce the systematic cosmic-shear errors below statistical errors for LSST 
r band. But <?-band dispersion effects remain large enough that it seems likely that in- 
duced systematics will dominate the statistical errors of both surveys, and cosmic-shear 
measurements should rely on the redder bands. 



1. Introduction 



Weak gravitational lensing (WL) can be used for high-precision measure ments of the expansion 
history of the Universe and the evolution of gravitational potentials within it ([Mellierlll999l ; iRefregier 
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20031 ; iHoekstra and Jainll2008l ) . The WL effect by large-scale structure is best detected as a coherent 
elongation of background galaxy images, "cosmic shear," which induces an RMS shift in axis ratios 
of only ~ 1% red or less for cosmologically distant source galaxies. Several ground-based surveys of 
> 1000 deg 2 are commencing this year, including the 5000 deg 2 Dark Energy Survey (DEsfl with 
the goal of measuring the amplitude and redshift dependence of cosmic shear to high precision. A 
yet larger and deeper survey project for late in this decade, the Large Synoptic Survey Telescope 
(LSST] 2 ], is under intensive development. 

Ground-based visible observations are subject to atmospheric refraction, shifting each photon 
toward the zenith by an angle R(X, z a ) that depends on wavelength and on apparent zenith angle 
z a . For finite-bandwidth observations, there are two measurable consequences: first, a centroid 
shift 

JdXR(X,z a )F(X)S(X) 

fdXF(X)S(X) [L) 

that depends on z a , the instrumental response function F(X), and the source spectrum 5(A). 
Additionally the dispersion stretches the image along the zenith, which can be quantified by the 
photon-weighted second central moment of the dispersion kernel: 

= JdX [R(X,z a )-R\ 2 F(X)S(X) 

JdXF(X)S(X) 1 j 

Like any other instrumental distortion or convolution of the galaxy images, this coherent vertical 
elongation of the galaxy images will be mistaken for WL shear if not removed. The standard strategy 
in WL observing is to use stellar images to estimate the instrumental point spread function (PSF), 
interpolating the stars' PSFs to the location of each galaxy, and then effect some limited form of 
deconvolution to extract the "pre-seeing" shape of the galaxy free of instrumental effects. In its 
simplest form, this instrumental correction amounts to subtracting the intensity-weighted second 
moments of the PSF from those observed for the galaxy to estimate the pre-seeing moments. If the 
atmospheric dispersion of the galaxies and stars are the same, then this stellar PSF correction will 
automatically correct the galaxy images for the dispersion Vo If) however, a galaxy's V differs by 
an amount AV from the mean stellar V used to estimate the PSF, then the inferred galaxy shape 
will be incorrect and propagate to cosmic shear errors. Our goal is to calculate the size of this 
effect, see if it will contribute a significant error to the cosmological results of DES or LSST, and 
evaluate some rudimentary techniques for eliminating the problem. 

Since cosmic-shear measurements are not affected by small shifts in the positions of galaxy 
images, we are less concerned with the centroid shift R of the galaxy image than with the broadening 



x http: / /www.darkenergysurvey.org/ 
2 http://www.lsst.org/lsst/ 

3 To be precise, the dispersion acts as a convolution, like the PSF, only if the galaxy's shape is constant across the 
band, i.e. the galaxy has homogeneous color. We will ignore in this work the complications of combining atmospheric 
dispersion with color gradients, as they produce shear errors at lower order. 
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V. Suppose, however, that we intend to combine information from exposures at multiple hour angles 
and/or filters, for instance by stacking the images before measuring shapes. The registration of 
images will normally be done by forcing overlap of stellar images. Spectral differences between the 
galaxy and the mean calibrating star will induce a difference AR in their mean refractions which 
depends on filter and hour angle. The stacked image will register the stars properly but not the 
galaxy, broadening the stacked galaxy image and hence producing a coherent shape error which 
again is confused with cosmic shear. Since both DES and LSST plan to combine many exposures 
to produce a single galaxy shape, we estimate AR and determine whether it causes cosmic-shear 
errors comparable to the statistical errors of these surveys. This problem can be avoided by leaving 
the galaxy centroid free to differ among exposures in a shape analysis — but this typically requires 
that the galaxy have high signal-to-noise ratio in a single exposure, which precludes use of most of 
the galaxies detectable with the full combined survey data. 

The effect of atmo spheric dispersion on weak lensing shear measurements was crudely estimated 



by lKaiser et al.l ()1995l ). who determined that it should be unimportant to the measurement of galaxy 
cluster masses. In the intervening 15 years, however, WL observers have become far more ambitious, 
aiming to measure the much smaller cosmic-shear signal to high precision, so it is necessary to re- 
evaluate these limits, especially sinc e neith er the DES nor LSST optics incorporate atmospheric 



dispersion correctors. iJee and Tvsonl (|201ll ) investigate the impact of atmospheric dispersion on the 
LSST observations, but only with the criterion that the dispersion V be a small contributor to the 
overall size of the PSF — they did not investigate the much stricter criterion that V be correctable 
to a small fraction of the expected cosmic-shear signal. 

In the next section we derive requirements on AR and AV^ such that they will not bias cosmic- 
shear measurements at the level of the statistical errors of the DES and LSST surveys. Then £j3] 
describes our methods for estimating these quantities, and gives the resultant values in the griz 
filters for a range of galaxy types and redshifts. Having demonstrated that atmospheric dispersion is 
indeed a significant issue in all but the z band, we ask in £}5] whether simple color-based corrections 
to AR and AV can recover the required accuracy. We conclude with an outlook in ^6l 



2. Tolerable levels of atmospheric dispersion 
2.1. Effect of dispersion on measured shear 



We use a simplified model of WL shear measurement to make rough estimates of the tolerable 
levels of AR and AV for cosmic-shear experiments. Galaxies are assigned shapes via their intensity- 
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weighted second central moments 



V = J J dxdyl(x,y)(fi - p)(v - u), (3) 

U = j J dxdyl(x,y)n, (4) 

f = J dxdyl(x,y) (5) 

where I(x, y) is the surface brightness distribution of the source. In practice these integrals have 
divergent noise so must be bounded by some window function, which complicates the measurement 
but has little impact on our estimates of the effect of the (small) atmospheric dispersion on the 
measured shapes. The ellipticities are assigned as 

ei - J^T (6) 

xx r yy 

and we can also assign the galaxy a second-moment-based radius, usually quite similar to its half- 
light radius, 

r 2 EE (I XX + Iyy) . (8) 

The two components (71,72) of the applied gravitational lensing shear are estimated as 27^ « (ei), 
where the average is taken over the ensemble of galaxies in a selected region of sky position and 
source redshift. 

These formulae require access to the galaxy's intensity I 9 (x,y) before any instrumental dis- 
tortions, but we only have access to the observed 1° after convolution with the PSF (including 
atmospheric dispersion components) I*. For present purposes it is a good approximation that: 

1. The second central moments of the observed image are the sum of those from the galaxy 
and the PSF, which means that the galaxy moments can be estimated as lj{ u = 1° — I*. 
This holds exactly for shape measurements with unweighted second moments, which are 
not practical. However the additivity of second moments under convolution also holds for 
Gaussian galaxies and PSFs, hence we can expect this relation to be a good approx i matio n 



for shape measurements based on Gaussian- weighted moments, such as lKaiser et al.l ( 19951 ). 



2. The effect of atmospheric dispersion is to boost slightly the second moment of the PSF in 
the zenith direction (call this the x direction), I* x — > I* x + V. This relation is exact for 
unweighted second moments, but is also a very good approximation for Gaussian PSFs in the 
limit when the atmospheric dispersion is much smaller than the PSF itself — which is certainly 
true here. The same boost is given to I° x , and the subtraction will cancel V in the estimate 
of intrinsic moments if our PSF is correctly calculated for the galaxy. 
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Under these conditions we see that a small error AV <C r 2 in the estimate of the V for a galaxy's 
PSF will result in an error in the zenith second moment, I% x — > 1% X + AV, propagating to shape 
errors 

/ AV\ AV 
ei ->• ei ( 1 - — J + — 



e2 



e 2 1 



r 
AV 



(9) 
(10) 



The shear estimate on a patch of sky is one-half the average galaxy ellipticity. In the shear- 
measurement literature it is typical to d istinguish a multiplic ative shear error m and an additive 
shear error c as 7? st = (1 + mi)7* rue + c, (jHevmans et al,ll2006l ). Comparing to the above, we infer 
that the atmospheric dispersion errors propagate to a multiplicative shear error 

'AV^ 



m 



(11) 



where the average is over the galaxies in a chosen bin of sky area and source redshift. The dispersion 
error also produces an additive shear error in the zenith direction of very similar amplitude, c\ ~ 
m/2. 



2.2. Shear measurement error budgets 



The multiplicative error m in shear measurement propagates directly to a fractional error 2m 
in the cosmic-shear power spectrum. Particularly pernicious are multiplicative shear errors that 
vary as a function of source redshift, since most of the power of cosmic-shear experiments comes 
from tomograp hic analyses of detai led comparison between shear strength on different bins of 
source redshift. iHuterer et al.1 ()2006l ) (hereafter HTBJ) calculate the degradation of cosmological 
constraints for the DES and LSST cosmic-shear surveys due to multiplicative shear calibration 
errors with unknown redshift dependence. From their Figure 4 we can infer that a systematic 
error in m of 0.008 (0.003) for DES (LSST) will raise the experiment's uncertainty in the dark- 
energy equation-of-state parameter w by a factor of v2 above the purely statistical errors. [It 
takes ~ 5x larger m values to double the cosmological errors above the statistical uncertainties.] 
Amara and Refregierl (j2008l ) carry out a similar analysis, concluding that |m| < 0.001 is required 
for a space-based cosmic-shear survey somewhat more ambitious than LSST, consistent with HTBJ 
at factor-of-2 level. Since atmospheric dispersion is just one of several potential sources of shear 
calibration error, an experimental error budget would likely need to hold the contribution from 
atmospheric dispersion to < i of the total value in HTBJ. We can hence set a requirement that 



< 



0.004 (DES) 
0.0015 (LSST) 



(12) 



Placing limits on additive shear errors is more complex because the limits on c depend strongly 
on how the additive signals correlate across redshift bins. We can make a crude estimate, however, 
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as follows: first, we note that the RMS value of cosmic shear is ~ 0.02, so the integrated shear 
power spectrum is ( r y 2 ) ~ 4 x 10 -4 . The maximum tolerable LSST multiplicative shear error of 
m < 0.003 produces an error in the integrated power spectrum of 2m{ / y 2 ) ~ 2.4 x 10 -6 . Spurious 
shear should not generate power larger than this value, (c 2 ) < 2.4 x 10 -6 , suggesting that we need 

crms < 0.0015. (13) 

If we again allocate half of this error budget to atmospheric dispersion, and keep in mind that the 
dispersion generates c\ = m/2, we find that the limits on AV^ induced by additive shear are very 
similar to those induced by multiplicative shear (Equation (|12p ). 



One issue that is more severe for additive shear than for multiplicative is that, for multi- 
plicative errors, we are on ly sensitive to the mean m over all galaxies in a given redshift slice 



(jGuzik and Bernstein! 120051 ). whereas for c it is the variance across the sky, not the mean value, 
that is of particular concern. So we should be alert to aspects of the atmospheric dispersion sys- 
tematic that will cause spurious shear fluctuating on the ~ 10' scales that carry most cosmic-shear 
information. For instance, if AV depends on source galaxy type within a redshift bin, we must be 
concerned that there will be variations in c across the field as the mean background type fluctuates. 



2.3. Requirements on atmospheric dispersion errors 

Our final step in setting requirements on AV is to estimate the galax y size to be used in 



Equation (fT2l) . The DES sets a goal of 10 lensing source galaxies per arcmin 2 . Ijouvel et al.1 (|2009l ) 
produce a model of the background galaxy population sizes and colors using HST/COSMOS and 
other observations; from this catalog we note that there are 10 galaxies arcmin -2 with half-light 
radii r > O'M, hence the largest possible value of (1/V 2 ) for DES would be to use this population. 
LSST is more ambitious, with a target source density of 30 arcmin -2 , which will require reaching 
down to galaxies with r > $.'27. For this latter population, (1/r 2 ) rj (O'.'SS) -2 . It is also true, 
however, that many of the galaxies above these cutoff r values will be unusable for weak lensing — 
too faint, crowded, or poor photo-z's — which will force the surveys to use galaxies smaller than 
these best-case cutoff sizes. As a rough guess, we therefore adopt estimates of (1/r 2 ) = (Cf.'4)~ 2 for 
DES and (0'.'27) -2 for LSST. The estimated requirements on AV for the target galaxy population 
are therefore 

' 6 x 10 -4 arcsec 2 (DES) 
1 x 10 -4 arcsec 2 (LSST) 



\^v)\<{ "";„_4 2 ioZ ( 14 ) 



Finally we consider requirements on AR, for galaxies where different exposures must be regis- 
tered for shape measurement. In this case the mis-registration will cause a blurring of the galaxy 
that will be roughly equivalent to an added dispersion AV ~ ((Ai?) 2 ). The details will depend 
upon the distribution of filters and zenith angles that one is trying to stack, but we can expect 
this variation to be of the same order as AR itself. So we set rough requirements that ((AR) 2 ) be 



-7- 



below the AV requirements in equation (|14|) . or: 

A - f 25 mas (DES) 

a-Rrms < { in / T ' 15 

10 mas (LSST) 

An additional point to make here is that two galaxies with opposite signs of AR are both elongated 
in the same direction after stacking, hence our requirement is not on the population mean of AR; 
it is on the RMS value of AR of individual galaxies. Galaxies with negative AR do not cancel the 
spurious shear of those with AR > 0. 

To summarize, atmospheric refraction has the potential to significantly bias the results of 
cosmic-shear measurements whenever the stellar calibration results in errors for galaxies exceed 
any of the following limits: 



1. The dispersion error (AV) of the galaxy population varies with redshift or other binning 
criteria by an amount exceeding 6(1) x 10 -4 arcsec 2 for DES (LSST). This produces a mul- 
tiplicative shear error that is too large. 

2. The galaxy-to-galaxy variation in AV within a redshift bin substantially exceeds the above 
values, which could produce an additive shear fluctuating on arcminute scales that biases the 
cosmic shear power spectrum. 

3. The galaxy-by-galaxy error in refraction centroid ARrms acquires typical value exceeding 25 
(10) mas in any redshift bin for DES (LSST). If this occurs, then stacks of images taken at 
different airmasses and/or filters will be misregistered and both multiplicative and additive 
shear errors will be produced. 



3. Atmospheric dispersion calculation methods 

To quantify the effects of refraction on astronomical images, we first write the refraction angle 

as 

R{X,z a ) = h(X)g(z a ) (16) 
where h(X ) and g(z a ) encode its dependence on wavelength and altitude angle. As indicated by 



Cox et al. 



(2001 



, g(z a ) can be approximated as g(z a ) = tan(z a ). This functional form implies that 
AR and AV scale as tan(z ) and tan 2 (z a ), respectively. We will present results assuming z a = 45°, 
since they are simple (tan 45° = 1) and the AV for a particular field of some survey should be scaled 
by (tan 2 z a ). The required scaling factor for DES is not, however, far from unity, so we will ignore 
it: we take a simulation of all the DES pointings produced by its scheduling software ObsTacq and 
bin them by declination. Figure Q] shows the relevant quantity (tan 2 z a ) vs declination expected for 



4 Provided by Eric Neilsen ( Neilsen 



2012h 
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DES, showing that (tan 2 z a ) is within about 20% of unity across the declination range of the survey. 
For LSST, Figure 3.3 of the LSST Science Boo^ suggests that the LSST i-band observations have 
(tan 2 z a ) ~ 0.6, so we can expect that the northern and southern extremes of the survey will have 
(tan 2 z a ) at or above 0.8, slightly below our nominal value of unity. 

The calculation of shear artifacts induced by AR in a stacking analysis is more complex, 
because AR is actually a vector quantity and the direction to zenith will vary with hour angle of 
observations of a given field. The deleterious effects will scale with the two-dimensional variance 
of the AR vectors over the course of the survey, so the relevant scaling factor will depend upon the 
hour-angle distribution as well as zenith-angle distribution, but will be < (tan 2 z a ). Our assumption 
of unity for this scaling factor is hence conservative, but we will see below that the shear systematics 
from AV are larger than those from AR even with this conservative upper bound on AR. 



We make use of values for h(X) tabulated in I Cox et al.l (|200ll ). slightly corrected to take into 
account the av erage conditions of pressure and temperature (770 millibar and 11 C, respectively 
(jStefanikl 120071 )) at Cerro Tololo, Chile (DES). As can be readily seen in Equations ([1]) and ([2]), 
R and V both depend on the spectrum of the source object, 5(A). To ma ke our calculations, we 
use the empi rical galaxy Spectra l Energy Distribution (SED) templates of IColeman et al.l (|l980l ) 
(CWW) and kinnev et al.1 (|l99d M The CCW templates consist of 4 galaxies of type E, Im, Scb, 
Scd and the Kinney templates are representative of both quiescent and starbust galaxies. For the 
star spectra, we utilize the stellar libraries provided by I Pickles! (|1998l ). We calculate the values 
of R and V for each stellar spectral type, and for galaxy types at values of redshift z E [0, 1] 
with Az = 0.02, according to Equations ([I]) and (J2]). Then, taking a G5V stellar template as our 
reference mean stellar spectrum, we calculate AR and AV for each object in each of the four griz 
filters and compare our results to the requirements for cosmic shear derived in £}2j 



4. Results for stellar and galactic spectra 

The top panels of Figures [2] and [3] show the calculated values of AR and AV in the g filter as a 
function of redshift (at z a = 45°). The remaining panels show these same quantities after applying 
a correction linear in color that will be described in the next section. Similar calculations and plots 
were performed and obtained in the r, i, and z filters. The results for stellar spectra at z = show 
how the effect of atmospheric refraction is greater for bluer sources, as expected. DES (LSST) 
requirements, as calculated in Equations (I14D and (I15p . are also shown in the plots. To determine 
whether these requirements are satisfied, we first examine the galaxies expected to be detected by 



5 http: // jhttp://www.lsst.org/lsst/scibook] 

6 These spectra are extended in the UV and IR regions by means of the synthetic Bruzual and Chariot 1 1993^ 
spectra. 
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the DES according to the DES Data Challenge 6B0 To each of the ss250,000 galaxies detected 
in a 10 deg 2 subsample of the simulation we assign a AV and AR from the template that most 
closely matches its colors and redshift. Then we divide the galaxy population into four redshift 
bins, computing the mean value and the dispersion, of AV and the RMS value of AR among the 
galaxy population. These are represented by the blue circular dots in Figures [6] and [7] (we assume 
that the LSST population has similar type distribution to the DES population). 

It can be seen that cosmic shear measurements performed in g band will, in the absence of 
some correction scheme, have systematic errors due to atmospheric dispersion that dominate the 
statistical errors, with AV up to 6 x our threshold of importance for DES and 30 x for LSST. The 
A-R rms values are also large enough to cause dominant systematic errors, but AV induces worse 
problems, so we will focus on AV henceforth. In r band, the levels of AV and A-Rrms remain 
large enough to dominate the LSST error budget but are 1.6x our nominal threshold for DES 
systematic errors, so may be acceptable for DES. For LSST, even in i band, the dispersion effects 
are up to 2.6 x larger than our derived requirements, so might be subdominant to statistical errors if 
our various approximations have been too conservative. But certainly the refraction effect requires 
attention in constructing an LSST error budget. Atmospheric dispersion is not a problem for either 
experiment in z-band (and we have verified that this is also true, as expected, in Y band). 



5. Strategies for calibrating effects of atmospheric dispersion 

The previous section shows that the effects of atmospheric dispersion on shape measurements 
are non-negligible in the g, r and i filters. As was pointed out above, AR and AV depend on 
the spectra of the observed objects, which cannot be measured in a photometric survey like DES. 
We would like to develop a calibration or correction that depends only on observable quantities, 
such as the color of the stars and galaxies. We calculate AR and AV of each star and galaxy at 
different redshift values and filters, and plot them as a function of color in Figures [H and [5l The 
plots suggest that a function of color can be defined that partially corrects the values of AR and 
AV. We determine AR and AV functions that depend linearly on color and minimize the RMS 
residuals of the stellar templates — this is an operation that could be done empirically from the 
observations. The stars in Figures [6] and [7J give the A-RrmS; an d mean and variance of AV, for the 
galaxy templates in 4 redshift bins after application of the star-based linear color correction to the 
galaxies. A linear function of color that minimizes the RMS residuals of AR and AV for the galaxy 
templates (crosses in Figures [5] and [6]) produces a calibration with smaller residuals for galaxies, 
but note that such a calibration would have to rely on theoretical calculations since AV cannot be 
measured from the data for the galaxies. The linear fits to both stars and galaxies are also shown 
in Figures 0] and \E\ The A-Rrms an d AV for the galaxies after applying the linear corrections at 
each filter are shown in Figures [6] and [7J Either linear correction brings the atmospheric dispersion 



Simulations created by the DES Data Management Project 



http:/ /cosmology.illinois.edu/DES / main / 
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Fig. 1 .— tan 2 z a for DES exposures binned by field declination. The exposure lists are produced 
by simulations of t he full PES s urvey using historical weather patterns and the ObsTac survey 
planning software ()Neilsenll2012l ). 
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Fig. 2. — The differential centroid shift AR as a function of redshi ft (Az = 0.02) in the g filter, 
for zenith angle 45°. Standar d templates of gal axies from CWW (jColeman et al.l Il980l ) (E, Sb, 
Sc, Im galaxies ) an d Kinney (IKinnev et al.lll996l ) (quiescent — Sa, Sb — and starburst — SB, SB6 — 
galaxies) were used. iPicklesI (119981 ) provides templates of stellar spectra ("ukxxyy" labels). A G5V 
spectrum was used as the reference to set AR = 0. The top panel shows the AR vaues for each 
spectrum at different values of z, whereas the lower panels show the same values after applying a 
linear correction in color, as described in the text. DES (LSST) requirements are that Ai?RMS fall 
within the light grey (dark grey) shaded region (Equation (fl5j) ). 
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Fig. 3. — The differential PSF second moment AV as a function of re dshift (Az = 0.02) in the 
g filter at zenith angle 45°. Standar d templates of galax ies from CWW (jColeman et al.lll980l ) (E, 
Sb, Sc, Im galaxies ) and Kinney (IKinnev et al,l Il996l ) (quiescent — Sa, Sb — and starburst — SB, 
SB6 — galaxies) were used. iPicklesI (|1998l ) provides templates of stellar spectra ("ukxxyy" labels). 
A G5V spectrum was used as the reference to define AV = 0. The top panel shows the AV^ vaues 
for each spectrum at different values of z, whereas the lower panels show the same values after 
applying a linear correction in color, as described in the text. DES (LSST) requirements are that 
AV fall within the light grey (dark grey) shaded region (Equation (I14p ). 
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Fig. 4. — AR as a function of color for stellar and galaxy spectra at different redshift values 
(Az = 0.02) for each of the griz filters at z a = 45°. The solid green (magenta) line represents a 
linear fit using only galaxy (stellar) spectra. Requirements for DES (LSST) on Ai?RMS are shown 
by the light grey (dark grey) shaded regions. 
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Fig. 5. — Ay as a function of color for stellar and galaxy spectra at different redshift values 
(Az = 0.02) for each of the griz filters at z a = 45°. The solid green (magenta) line represents a 
linear fit using only galaxy (stellar) spectra. Requirements for DES (LSST) on AV are shown by 
the light grey (dark grey) shaded regions. 
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effects within requirements in the r (i) band for DES (LSST). The LSST r band systematic errors 
are still above our thresholds of concern after linear color corrections, by factors up to 5 depending 
on the redshift bin and type of correction — thus must still be considered an important contributor 
to the LSST error budget, although suggestive that more sophisticated correction schemes could 
bring atmopheric-dispersion systematics back below statistical errors in this case. For g band, 
however, the linear corrections still leave errors well above our thresholds in both DES and LSST, 
and in fact the "corrected" data are worse in most of the redshift bins than uncorrected ones. Also 
we note that the galaxy-to-galaxy standard deviation of AV within a redshift bin is typically much 
larger than the DES or LSST requirements we have derived, suggesting that spatial variations of 
the induced additive shear will be an issue. It appears that atmospheric dispersion is likely to 
create dominant systematic errors in g band, and the surveys should plan on basing cosmic-shear 
measurements on the redder bands. We have verified that u band is, as expected, even worse than 
9- 
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Fig. 6. — Ai?RMS i n f° ur equally-spaced redshift bins for each of the four griz filters when no 
correction (blue dots), a linear correction in color using only star spectra (magenta stars) and a 
linear correction using only galaxy spectra (green crosses) is applied. DES (LSST) requirements 
are that Ai?RMS fell within the light grey (dark grey) shaded region (Equation (|15j) ) . 
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Fig. 7. — AV in four equally-spaced redshift bins for each of the four griz filters when no correction 
(blue dots), a linear correction in color using only star spectra (magenta stars), and a linear cor- 
rection using only galaxy spectra (green crosses) is applied. The vertical bars depict the standard 
deviation of AV between different galaxy templates and redshifts within each bin. They therefore 
represent the range of variation that the mean A V might accrue due to variations in the population 
mix of galaxies with redshift and with position on the sky. DES (LSST) requirements are that AV 
fall within the light grey (dark grey) shaded regions (Equation (114ft ), 



6. Conclusion 



Atmospheric dispersion can have a significant effect in PSF correction of shape measurements 
in broadband photometric surveys by affecting the first and second moments of the observed im- 
ages. If uncorrected, it produces cosmic-shear spurious signals larger than the desired statistical 
uncertainties for the g and r bands of the DES, and for the g, r, and % bands for LSST. This 
effect depends on the spectrum of the source and therefore is not easy to correct. A simple linear 
correction to the PSF position and size based on the observed color, calibrated empirically from 
stars or estimated theoretically for galaxies, should improve calibration of dispersion effects signifi- 
cantly. Even after applying this correction, however, weak lensing measurements in the g band are 
significantly degraded by dispersion effects for both DES and LSST. Cosmic shear measurements 
can be performed with images taken in the r band for the accuracy required by DES, but further 
work will be needed in the case of LSST. Shape and shear measurements in the z and Y band can 
safely be performed without attention to dispersion corrections. 

There are several possible tactics for further reduction of cosmic-shear systematic errors from 
atmospheric dispersion: 



• Use an atmospheric dispersion corrector (ADC) during the imaging. While this option is 
available for KIDS cosmic-shear surveys being conducted on the VLT Survey Telescopd§ and 
for the HyperSuprime Camera to be installed on SubarvJ§, an ADC is not part of the optical 
designs for DES and LSST. 

• Observe at lower zenith angles — AR and AV scale as tan(z a ) and tan 2 (z a ), respectively, so 
would be reduced by factors of y/3 and 3, respectively, if we were to take nominal zenith 
angles of 30° instead of the 45° assumed in our calculations. In practice it would not be easy 
to maintain such low airmasses throughout a large survey, especially if one wants to survey 
15,000 deg 2 or more of low-Galactic-latitude sky from the latitudes of ~ —30° of both Cerro 
Tololo (DES) and Cerro Pachon (LSST). 

• Derive a more elaborate correction than the linear fit we have used, e.g., an approach that 
makes use of multi-color photometry to estimate the full spectrum of each galaxy — which 
of course is already done in the course of most template-fitting estimators for photometric 



redshifts. This technique is discussed by ICvpriano et al.l (|2007l ) in the context of calibrat 



ing wavelength-dependent PSF variation across the very broad band planned for the Euclid 
spacecraft. In the Euclid case it is telescope diffraction, not atmospheric refraction, that 
causes the wavelength dependence of the PSF, but the issue is similar. 

One could simulate the full effect in detail to estimate additive and multiplicative corrections 



http://www.eso.org/public/teles-instr/surveytelescopes/vst.html 
9 http: / /www.naoj .org/ 
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that must be applied to the measured shear fields. This approach would have the usual 
limitations of precision calibration by simulation: it would be accurate only to the extent 
that the simulation captures the full multivariate distribution of the sizes, shapes, redshifts, 
and SEDs of galaxies in the true sky. The SED distribution will be particularly challenging 
to constrain for the faint end of the surveys' target galaxy populations. 

It seems likely, therefore, that both DES and LSST can make use of r or redder bands for cosmic- 
shear surveys, but this will require explicit galaxy-by-galaxy correction for PSF and astrometric 
dependence on the source spectrum. For DES, we show that a simple linear color correction will 
suffice, but more work is needed to devise a correction that recovers the full cosmic-shear utility of 
the LSST r band data. 

We note lateral color in telescope optics causes spurious spectrum-dependent cosmic-shear 
signals in a manner very similar to the atmospheric dispersion that we examine here. As with 
atmospheric dispersion, it is typical to specify an optical design that keeps this contribution well 
below the size of the PSF, but this is not a sufficient condition to guarantee that differences between 
galaxy and stellar spectra produce PSF calibration errors below the requirements of cosmic-shear 
surveys. 

This work was supported by NSF grant AST-090827 and DOE grant DE-FG02-95ER40893. 
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